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Abstract 

Adopting the thin-layer improved brick-wall method, we investigate the thermodynamics 
of a black hole embedded in a spatially fiat Friedmann-Lemaitre-Robertson- Walker universe. 
We calculate the temperature and the entropy at every apparent horizon for arbitrary solution 
of the scale factor. We show that the temperature and entropy display a non-trivial behavior 
as a functions of time. In the case of black holes immersed in universe driven by phantom 
energy, we show that for specific ranges of the equation-of-state parameter and apparent 
horizons the entropy is compatible with the D-bound conjecture, even the null, dominant and 
strong energy conditions are violated. In the case of accretion of phantom energy onto black 
hole with small Hawking-Hayward quasi-local mass, we obtain an equation-of-state parameter 
in the range w < —5/3, guaranteeing the validity of the generalized second law. 
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1 Introduction 



The discovery that the current universe is in an accelerating expansion phase, obtained from 
the observations of type la supernovae [U [2], inaugurate an era of intense theoretical research 
to understand the mechanism driving this accelerating expansion. An immediate consequence is 
that gravity behaves differently on cosmological distance scales. A variety of possible solutions to 
the cosmic acceleration puzzle have been debated during this decade including the cosmological 
constant, exotic matter and energy, modified gravity, anthropic arguments, etc. The most favored 
ones are dark energy based models and modified gravity theories such f(R) gravity and DGP 
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gravity with an equation of state parameter w = Pj g < — 1 (where p and P are the energy density 
and pressure of the cosmic fluid, respectively). 

The dark energy component is usually described by an equation-of-state parameter w < —1/3. 
The simplest and traditional explanation for dark energy is a cosmo logical constant, which is usu- 
ally interpreted as the vacuum energy. However, the value required to explain the cosmic expansion, 
which is of the order of 10~ 120 , cannot be explained by current particle physics. This is known as 
the cosmological constant problem. On the other hand, the first year WMAP data combined with 
the 2dF galaxy survey and the supernova la data favor the phantom energy equation-of-state of 
the cosmic fluid w < — 1 over the cosmological constant and the quintessence field. A candidate 
for phantom energy is usually a scalar field with the wrong sign for kinetic energy term [3], 0]. 
One crucial fate of an expanding universe driven by phantom energy is that the phantom energy 
density and the scale factor diverge in finite time, ripping apart all bound systems of the universe 
(galaxies, stars, atoms, nuclei), before the universe approaches the Big- Rip singularity [3J [5] . 

In this work we are interested by the effect of cosmological expansion on local systems, namely a 
black hole embedded in a an expanding Friedmann-Lemaitre-Robertson- Walker (FLRW) universe 
driven by phantom energy. The study of non-stationary event horizons and their thermodynamical 
parameters, relevant for a quantum theory of gravity, are currently attracting a great interest [6]- 
[19J. In the pure expanding flat FLRW universe dominated by phantom energy, the radius of 
the observer's event horizon diminishes with time and consequently the horizon entropy, Sh < 0. 
An other example is the quasi-de Sitter space, where the event horizon and apparent horizon are 
different [201 [2Tj . It was found that the first law and second law of thermodynamics cannot hold 
when both the horizons are considered. Other authors discussed the effect of the presence of black 
holes in expanding phantom energy-dominated universe on the validity of the generalized second 
law of gravitational thermodynamics (GSL) [221 123] . 

In this paper, we present a detailed calculation of the entropy of the solution of Einstein's 
equations recently found and representing a black hole embedded in an expanding Friedmann- 
Lemaitre-Robertson- Walker (FLRW) universe j24l [25] . Our approach is based on the brick-wall 
method (BMW) invented by t'Hooft [26J. In this method a brick wall around the horizon prevents 
the divergence of the free energy (or entropy) of the black hole, which is identified with the 
canonical ensemble statistical-mechanical free energy (or entropy) due to quantum excitations in 
thermal equilibrium with the black hole. However, the BMW can not be used for non-equilibrium 
systems, like black holes with multi-horizons, where each horizon can be considered as an isolated 
thermodynamical system. The solution to the problem, known as the thin-layer improved method, 
consists in considering a thin layer near the horizon as a local equilibrium system and invoking 
thermodynamics of composite systems, such that the entropy of a multi-horizons black hole is the 
sum of the contributions arising from each horizon [27]. Then in the improved BWM, the global 
equilibrium has been replaced by a local equilibrium on microscopic scales. This method has been 
applied to a wide range of black holes with multi-horizons like Schwarzschild-de Sitter black hole, 



Kerr-de Sitter black hole, Vaidya black hole, the 5D Ricci-flat black string and recently to the 
calculation of the entropy in the pure expanding FLRW universe |28| . 

The organization of the paper is as follows: in Section II we review the exact solution describing 
a black hole embedded in an expanding Friedmann-Lemaitre- Robertson- Walker (FLRW) universe. 
In section III, the statistical mechanical free energy, thermal energy and entropy due to a scalar field 
are computed using the improved thin-layer BWM. Then we discuss the conditions under which the 
D-bound conjecture [29J is protect in expanding universe in the presence of black holes. In section 
IV, our attention will be focused on the conditions under which the second law of thermodynamics 
and the GSL are satisfied when the universe is driven by phantom energy. Finally, we discuss and 
summarize our results in section V. 



2 Cosmological expanding black hole 

The first solution of Einstein's theory of general relativity describing a black hole like object 
embedded in an expanding universe was introduced by McVittie in 1933 [30j, and is given in 
isotropic coordinates by 

( l Mq \ 2 4 
ds 2 = _\ M<K) dt 2 + a 2 {t) ( 1 + _M0 \ ( dr 2 + r 2 d ^ (1) 

\ l + 2a(t)r ) 

where a(t) is the scale factor and Mq is the mass of the black hole in the static case. In fact, when 
a(t) = 1, it reduces to the Schwarzschild solution. When the mass parameter is zero, the McVittie 
reduces to a spatially flat FLRW solution with the scale factor a(t). The global structure of (CQ) has 
been studied and particularly it has been shown that the solution possesses a spacelike singularity 
on the 2-sphere r = M /2, and cannot describe an embedded black hole in an expanding spatially 
flat FLRW universe [3HE2]. 

In the following we adopt the new solution describing a black hole embedded in a spatially flat 
FLRW universe [21 [25] 

ds 2 - j^l^^ 2 + (l + ^f) + ^) > ( 2 ) 



Using the areal radius 
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the metric takes the following suitable Painlevele-Gullstrand form 
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dtdR + R 2 dn 2 , 



where H = a/a is the Hubble parameter and overdot stands for derivative with respect to the 
cosmic time. We observe that the term R 2 H 2 plays the role of variable cosmological constant. In 
order to write the metric in the Nolan gauge, we introduce the time transformation t — > i to 
remove the dtdR term 



dt = F~ l (t, R) 



dt + 
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where the integrating factor F (t, R) satisfy 
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Then substituting (j3J) into fflh and replacing t — > t , we obtain 
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ds A = dh z + R 2 dn 2 



where the two-dimensional metric is 
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The apparent horizons (AH) are solutions of the equation h ab d a RdbR = 0, which leads to 

2m H (t) 



/? T RH)\ Ra =0, 
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dR 2 . (8) 
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where we have used the Hawking-Hayward quasi-local mass m^(t) = M a(t). This mass is always 
increasing in an expanding universe. Discarding the unphysical branch with the lower sign, we 
obtain 



R 



c 



2H 



lrh H (t) 



Rb ~ (l - y/1 ~ 8m H (t)) , 



(10) 
(11) 



where Rq and Rb are the cosmic and the black hole AH, respectively. We observe that the AH 
coincides at a time t* for which d(i*) = 1/8M . This coincidence takes place in a future or past 
universe depending on the kind of matter accretion onto the black hole. 



3 Statistical mechanical entropy 

Now we consider the statistical mechanical entropy that arises from a minimally coupled quantum 
scalar field of mass \x in thermal equilibrium at temperature /3 using the thin layer improved BWM 
method. 

The field equation on the background given by (BJ is 

(□ - /, 2 ) $ (x») = 0, (12) 

where $ (a;") = (R, t) Y lm (9, <p). 

Using the relation 2m# (t) = Ra (1 — RaH), the Klein-Gordon equation takes the form 
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V(R,t) = 0. 



In the semi-classical approximation, the wave-function is obtained by the ansatz \1/ (R, t) ~ e l ^« h n s n (R,t) _ 
Neglecting higher orders contributions and substituting in Eq. (1131) with the ansatz Sq = —iut + 
iB(R,t), we obtain 
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In our derivation we have assumed the constancy of the frequency u near the AH and that 
B(R,t) <C u. Solving for B' we obtain 



b'(r, o s Pi- - m (1 _ _ RaH)) ± j^^m - m w + *m + d). us) 

where we have set 



f(R) = i - iM (i _ 
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which is just g tt . The sign ambiguity in Eq. (|15|) is related to outgoing (dS^/dR > 0) or ingoing 
(dS /dR < 0) particles. 

The number of radial modes with energy less that u is defined by 
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UJ ' l ) = WZ / dR (pt ~ Pr) 



(17) 



where we used the average of the radial momentum. The second term in (TTT1) is caused by a 
different direction such that the outgoing and ingoing particles are taken into account. Using 
Eq. ffT5]) we get 
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Now, according to the canonical ensemble theory the free energy is defined by 



(3F = J dN\n(l-e-^) 



Integrating by parts we obtain 
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where N (u) is the total number of modes with energy less than u given by 
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the free energy takes the form 
F = - 
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N(w) = j dl(2l + l)n(u,l). (21) 
Substituting the expression of n and restricting the integration over / to the range 
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The integral over the radial variable is determined by the improved thin-layer BWM boundary 
conditions 



$ (x») = for R B + h B <R<R B + h B + 6 B , (24) 
$ = for R c - h c - 5 C < R < R c ~ h c , (25) 

where h B ,hc <C R B , Rc are ultraviolet cutoffs and S B , 5c are thickness of the thin layer, near the 
BH horizon and cosmomlogical horizon, respectively. Then, the integrals over R are performed by 
noting that f(R) — > in the near vicinity of the AH. In fact f(R) can be approximated by 

/ « -|- (1 - 2R A H) (R - R A ) . (26) 
Then using the near horizon approximation (u 2 — fi 2 f(R)) 3 ^ 2 — > u 3 , and the following integral 

r^-=^- (27) 

the free energy is written as F = F B + Fc with 
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where (3 B and /?c iare the Hawking temperatures evaluated at the BH and cosmic AH, respectively. 
Using f[2"61) we obtain 
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Let us now introduce the invariant distances h B and 5 B defined by 
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The invariant distances near the cosmic AH, he and 5c, are obtained by the substitutions (h B , 5 B ) — 
(he, 5c) and (1 — 2R B H) — > {2R C H — 1) in the above relations. Then, the free energy takes the 
form 
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The internal energy and entropy defined by U 
then obtained as 



A (/3-F) and S = f3 (U — F), respectively, are 
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Using the area of the BH horizon and cosmic AH, Ab = 4:irR B and Ac = 4:7rR c , respectively, we 
rewrite the entropy as 
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Now the inverse temperature is calculated at the AH using /3 1 = \ka \ /27T, where the dynamical 
surface gravity associated with the AH is ka = (2y/^K) d a (y/^hh ab d b R) \ Ra and where the 
metric h a b is defined by ds 2 = h ao dx a dx b + R(x)dQ 2 . The calculation yields 
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where H = a(t)/a(t), H = dH/dt. We also show that the Misner-Sharp energy defined by 



M {R, t) = f (1 - g^dpRduR) is given by 
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Writing nin is terms of the AH, the temperature and the Misner-Sharp energy become 

1 
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The temperature associated with dynamical horizons is not only related to the Misner-Sharp mass 
on the AH, but also to the matter content of the universe [33]. Here, we point that in a series of 
papers devoted to thermodynamics of expanding FLRW universe with dark energy it is assumed the 
usual law T4 = 1/2ttRa for stationary horizons[34J, which certainly invalids the results obtained 



[35J. The expression of temperature (1391) is similar with the temperature in York theory of black 
holes confined in isothermal cavities [36J. In this context, a local observer at rest will measure a 
local temperature which scale as 1/ y/—gu f° r a self-gravitating mass in thermal equilibrium with 
the boundary of the cavity. The effect of the boundary is reflected in the relation (13811 . which state 
that the Hawking-Hayward quasi-local mass is the sum of the thermal energy and a gravitational 
self-energy due to the cavity. In our case, the black hole is embedded in an expanding universe, and 
the cosmic apparent horizon plays the role of the boundary of the cavity. On the other hand the 
behavior of temperature with time is contrasted with the naive expectation that the temperature 
is redshifted by the scale factor a(t) relative to its value in static background. The additional 
correction terms to the temperature become important when the density energy of the universe 
is negligible compared to the expansion rate. The same remark apply to the entropy which will 
show a behavior with time different from the expected one S ~ S 'static® -3 ■ Here some comments 
are appropriate about the assumption of local equilibrium used in the calculation, even we are 
considering dynamical horizons. The improved thin-layer BMW method is based on the fact that 
Hawking radiation is mainly due to the vacuum fluctuations near the black hole horizon(s), and 
then allows to identify the black hole entropy with the statistical-mechanical entropy associated 
with quantum fields outside the black hole and confined in small region(s) in the vicinity of the 
horizon(s). This cure the problem of maintaining thermal equilibrium at large scales of the original 
BWM of t'Hooft, without taking into account the back-reaction of the background. However, in 
small regions we are dealing with the assumption of local equilibrium and the validity of the 
statistical laws we used in the calculation, and therefore need a physical justification. The local 
equilibrium state is achieved if thermodynamics properties of the system are slowly varying. The 
last condition when applied to the case of the case of a black hole immersed in an expanding 
FLRW universe is equivalent to the following statement: "local equilibrium is achieved if local 
thermalization of the system is faster than the expansion rate of the universe". In fact assuming 
that if <C 1, and using Eq.(39) we obtain 4jr ~ ^p- — ^jf <C 1. This condition can be rewritten as 
jj ^> 1, which is an expression of the last statement 

In the pure flat FLRW case where only the cosmic AH survives, Rc = 1/if, we obtain 
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which is proportional to the de Sitter temperature, [3 c l = H/2tt. 
Now, substituting Eq.( !39|) in (|36|) , the entropy takes the form 
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In order to compare our results with the standard FLRW universe, we set the cutoffs as 



and introduce the following factors 
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Then the final expressions of thermal energy and entropy take the forms 
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We assume now that matter in the universe is in the form of an imperfect fluid with a radial 
heat flux, described by the stress-energy tensor 
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where = 0, 0, 0) is the fluid four velocity, = (0, q, 0, 0) a spacial vector describing radial 
heat current, and A = (l + M /2r), B = (1 — M /2r). The general solution of the Einstein 
equations are given by 
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Assuming r ^> m and radial energy inflow {q < 0), the accretion rate can be written as 

mn — GaA \q\ , 

where A = J J dOdip^/g^ = 4TTr 2 a 2 A A . In terms of the comoving AH, these equations reduce to 
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At this stage we would like to comment about the law of temperature T4 = 2 J- R , widely used in 



the literature. With the aid of Eqs. (15011511) we observe that the later law is only valid if the energy 
density of the universe is significant relative to the expansion rate 



p > H. (52) 

Now, it is instructive to rewrite the dynamical surface gravity as an invariant. In fact using 
the Misner-Sharp mass M (Ra,T) = Ra/2 and defining the projection of the (3+l)-dimensional 
stress-energy tensor in the normal direction of the 2-sphere, T A ^ = h ab T a i, (a, b = 1,2), we obtain 
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where we have set G — 1. The later invariant relation has been recently proved in different gauges 
|37| . On the other hand this equation leads immediately to the analogue of first law for dynamical 
black holes. In fact using A a = ^nR\ and Va = \kR? a as the AH area and volume associated with 
the AH, we obtain 

dM = ^dA A - \r A 2) dV A . (54) 

Assuming now an equation of state of the form p = wp and using Eqs. (j50|5ip . we rewrite the 
factors J~b c as functions of the EoS parameter w, 
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where we defined the reduced AH Xb,c — Rb,cH. The factor in the standard FLRW universe, 
already given in (j4"Tl) . follows by setting xq = 1 in the second relation 



l-3w 
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which is exactly the entropy factor obtained in [28]. In this case the temperature associated with 
the cosmic AH becomes 

Pc 1 = (57) 

Since Tq > 1 in a universe driven by phantom energy, this result put some doubts on the claim 
that the GSL in the presence of a black hole in a phantom energy-dominated universe is protected 
if the temperature is written as f3^ 1 = bH/2ir, with < b < 1 [38]. In this work, the backreaction 
of the phantom fluid on the black hole has been ignored as the authors assumed the usual pure 
flat FLRW metric. 

In figure 1 we plotted the behavior of the factors J~b and J~c as functions of the EoS parameter 
w for different values of the BH and cosmic AH, respectively. It is interesting to note that the 



factors Tb and Tc do not depend on the AH when w = 



-1/3. 
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Figure 1: Variation of the entropy factors Tb and Tc with the EoS parameter w for different 
reduced AH Xb and Xc on the left and right panel, respectively. 



In the following let us check if the expression of the total entropy meets the D-bound conjecture 
|29j generalized to dynamical black holes. The later conjecture asserts that the entropy of a system 
within a boundary is less than or equal to the gravitational entropy, S < Sc = A/ AG. We consider 
the special ranges of values of Xb and Xc used in figure 1 , and look for the possible values of the 
EoS parameter w for which the D-bound conjecture is protected. In fact, from the relations giving 
the AH in Eqs. (llOHTTj) . and assuming that Mq ^ 0, we have < Tc < 1 for 1 < xc < 0.6 and 
-0.55 < w < 1/3, and < T B < 1 for < x B < 0.4 and -4/3 < w < 0. Then, within these 
ranges of the EoS parameter and the AH, it is immediate to verify that the entropy associated 
with the BH and cosmic AH satisfy the D-bound conjecture, 
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respectively. When Mo = 0, the BH horizon disappears and we obtain Sc < ^ for — 1 < w < 1/3. 
To verify the D-bound conjecture for the total entropy Sb + Sc, we have to verify that the following 
condition 



T B R% + T C R 2 C < Rl (59) 
holds, where Rq = 1/H corresponds to the AH in the pure flat FLRW universe. Using i?o = 



Rb + Rc, the condition (1591) is rewritten as 



J- c <l + (l-^)4 + 2— . « 60 ' 
where we have used the reduced AH. Since we have Xb/%c < 1, the condition (160]) becomes 

Tc + ?b< 4- (61) 

We easily verify that the above condition is verified for an EoS parameter in the interval —0.8 < 
w < 0.6, and < x^ < 0.4. and 1 < xq % 0.6. On the other hand, if we restrict the values of the 
reduced cosmic AH to 1 < xc < 0.8, the condition (IBTj) is satisfied for —1.35 < w < 1.7. The later 
values of xc corresponds to a small Hawking-Hayward quasi local mass. Finally, we point that even 
the phantom field is unstable and violate the null, strong and dominant energy conditions (Null, 
SEC, DEC) [39j, the entropy of dynamical black holes immersed in phantom energy-dominated 
universe, meets the D-bound conjecture. 

In what follows, we would like to study the time evolution of the thermodynamics parameters 
when the black hole is embedded in an FLRW universe driven by phantom energy and accreting 
this cosmic fluid. In this situation, the scale factor is given by 

a ph (t) = a (t np -t)^+v , (62) 

where t r i p in (1B"2"|) is the big rip time. In that time increases, the black hole horizon 

increases monotonically while the cosmic AH decreases monotonically. The AH coincide when 

| 3(m + l) | 

, and after which the black hole singularity will become naked in a 
finite time, violating the cosmic censorship conjecture [25]. In figure 2 we show the evolution of 
entropy and free energy with time for w = —1.2, and different initial mass of the black hole. We 
observe, that entropy is decreasing and increasing at early times and late times when approaching 
t*, respectively, and diverges at the critical AH radius Rc = 1/2H. This indicates that at times 
approaching if:*, the second law of thermodynamics is protected, S > 0, while at early times the 
second lwa of thermodynamics is violated. The important result to note here is that, taking into 
account the effect of backreaction of the phantom fluid on the black hole in an expanding universe, 
the second law of thermodynamics becomes valid when approaching the coincidence time, while it 
is always violated in a purely phantom-energy dominated universe. 
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Figure 2: Variation of entropy and free energy as functions of time with w = —1.2. From right to 
left we have M = (solid) , 0.01 (dot) , 0.1 (dash) , 1 (dash-dot) . 



4 Generalized second law 

It has been advanced that when a black hole is embedded in an expanding universe driven by 
phantom energy, the GSL is verified under some restrictive conditions [22| 123] . Let us now proceed 
to discuss the GSL, S + Sf > 0, where S is the geometric entropy associated with the AH and Sf 
the entropy of the cosmic fluid confined between the BH and the cosmic AH. The entropy of the 
fluid can be obtained by using the Gibbs equation |40] . 



where T is the temperature of the energy in the vicinity of the AH. Using the continuity equation, 
which is given in our model by 



and assuming that the cosmic fluid is in thermal equilibrium with the boundary, the evolution of 
the fluid entropy is then given by 



TdS fA = Vdp +(p + p) dV, 
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Using this relation and since the temperature is meaningful only in the near horizon region, we 
assume that the fluid entropy is the sum of the contributions near the BH and cosmic AH, respec- 
tively, 
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A=B,C 

Let us now consider explicitly the case where the Hawking-Hayward mass is enough small so 
that the quantities associated with the BH horizon can be neglected. Then, we can set Rb — 
and the horizon entropy is due essentially to the contribution near the cosmic AH 



Sc ~ J^qtcRq, 

where J-q is given by (15611 . The time derivative of the phantom fluid entropy becomes 
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Now taking the time derivative of Sc, we obtain 



Sc + Sf 



R 2 C H 



2T ( 



c 



RcH 
RcH 



H 

1 + 2 jp 



+ 2F nR c R. 



c- 



(67) 



(68) 



(69) 



The GSL states that Sc + Sf > 0. Note that in the phantom era H > 0, Rc < 0, and since Tc > 0, 
a necessary condition for the GSL to be satisfied is 
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RcH 
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1 + 2 



H 

W 2 



> 0. 



(70) 



This condition can be integrated and leads to 



Rc{t) < a(£)2[w^]. 



(71) 



Using the expression of the cosmic AH, we get the following condition on the derivative of the 
Hawking-Hayward quasi-local mass 
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(72) 



Now substituting the scale factor given by (162]) . and the expression of the Hubble parameter we 
finally obtain 



1 

m H > m H ,crit = o 

o 

In figure 3, we plotted the variation of rriH,crit with time for different values of the EoS parameter. 
Knowing that the Hawking-Hayward mass is an increasing function of time in an expanding uni- 
verse, we observe that the GSL remains protected from the past to the present time if w < —5/3. 
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Figure 3: Variation of the derivative of the Hawking-Hayward quasi-local mass as a function of 
time. 
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Using rfiH = M a(t) and considering the present time, we get for the mass of the black hole 
the following constraint 

„, ^ 1 5 + 3w 
M > k. 

Assuming the positivity of the mass, we obtain again w < —5/3, in order for the GSL to be satisfied. 
This value of the EoS parameter is compatible with the analysis performed on the validity of the 
D-bound conjecture in the case where the cosmic AH is close to the critical value 1/2H. 



5 Conclusion 

In this paper we have studied the thermodynamical properties of black holes immersed in an 
expanding spatially flat FLRW universe. We have particularly calculated the entropy and temper- 
ature associated with the apparent horizons using the improved thin-layer brick wall method and 



the dynamical surface gravity, respectively. The temperature and entropy at the apparent hori- 
zons (AH) display a non trivial behavior as a function of time, and are not scaled by the expected 
factors a(t) and a~ 3 (t), respectively. The correction terms become relevant when the expansion 
rate is significant relative to the density energy of the universe. On the other hand, we found that 
the sum of entropies associated with the AH meets the D-bound conjecture for an EoS parameter 
in the interval, —1.35 < w < 1.7, although for w < —1, the null, strong and dominant energy 
conditions are violated. We have also discussed the validity of the second law of thermodynamics 
and the generalized second law for a black hole embedded in phantom energy-dominated FLRW 
universe. The analysis showed that for arbitrary Hawking-Hayward quasi-local mass, the second 
law of thermodynamics is protected when approaching the coincidence time, at which the appar- 
ent horizons coincide, Rb = Rc — 1/2H. On the other, in the case of small Hawking-Hayward 
quasi-local mass, the GSL is only satisfied if w < —5/3. 
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